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16, Abstract 

This report presents the results of a modal transient analysis of the 
engine/- aircraft system. The response of the JT9D to analytically simulated 
vertical gusts arid landings was predicted using a NASTRAN (NASA STRUCTURAL 
ANALYSIS) finite element mathematical model of the JT9D/747 propulsion system. 

The NASTRAN finite element model of the propulsion system included engine 
structural models of the fan, 1 ow-/high-pressure compressors, diffuser/tur- 
bine cases, and high-/low-pressure rotors, as well as nacelle models of the 
inlet cowl, tailcone, and wing pylon. 

The analysis conducted predicts that an insignificant level of JT9D-7 perfor- 
mance deterioration would occur due to a typical vertical gust encounter or a 
typical revenue service landing. Analysis of a high sink rate landing with a 
heavy fuel load indicates the possibility of local wear, however, the lack of 
an accurate dynamic rotor/seal interference model precludes an accurate quan- 
titative evaluation of performance change for this once-per-airframe-1 if e 
event . 
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SECTION 1.0 
SUMMARV 


A jet engine's efficiency Is strongly dependent upon clearances between 
rotating blades and the corresponding stationary gas-path sea s. Ana- 
lytical studies examining the effects of various flight cycle loads 
upon powerplant clearances Suggest that rotor/seal 
tribute significantly to the short-term performance defloration ®J" 
hlblted by high bypass propulsion systems. In prior analytical studies 
m* ** flight cycle loads were treated In a quasi* steady, static manner, 
even though some flight cycle loadings are decidedly dynamic In nature. 
Typical examples of dynamic flight loads are gust 
landings. To better quantify the structure response ° f a %r*S B rS D Sl- 
f light environment, a modal transient analysis of the alrframe/propul 
sion system has been performed. The response of the JT9D-7 to analyti- 
cally simulated vertical gusts and landing touchdowns was predicted 
using an existing NASTRAN 9 (NASA STRUCTURAL AN alysis ) fin1tg glenient 
mathematical model of the JT907747 propulsion system. 

The NASTRAN finite element model of the propulsion system included en- 
gine structural models of the fan, low-/high-pressure compressors, dif- 
fuser/turbine cases, and the h 1 gh-/low-pressure rotors, as well as 
nacelle models of the inlet cowl, tailcone, and wing pylon. A modal 
synthesis technique was used to extract the integrated propul slon/air- 
frame system natural frequencies. Calculated frequencies were in flood 
agreement with available test data. Transient response cal cul at o 
included rotor gyroscopic terms and modal damping. Time dependent 
powerplant displacement data was related to performance characteristics 
usihd a postprocessing system allowing clearance changes to be 
ally reviewed and recorded at various time points during the load tran- 
sients. 

The analysis conducted predicts that an insignificant level °f /^D~7 
performance deterioration would occur due to a once-per-flight ver- 
tical gust encounter (0.04 percent ATSR at cruise) or a typical reve- 
nue service landing (0.2 percent ATSFC at cruise). These results 
strengthen the conclusions from the previous study (1) that steady 
aerodynamic pressure loads are probably responsible for the bulk of the 
TSFC loss predicted to be caused by flight loads. The results also 
indicate that the quasi-steady state approach to flight loads modeling 
is adequate to investigate the factors important to the deterioration 


* Numbers in ( ) indicate references in Section 9.0. 

** The once-per-flight gust encounter simulates response of the air- 
craft to a column of rising air and represents the maximum effects 
of air turbulence expected to occur during an average flight. 


process, especially the short^rm aspect. Analysis of a high sink rate 
landing with a heavy fuel load Indicates the possibility of local wear, 
however, the lack of an accurate dynamic rotor/seal Interference model 
precludes an accurate quantitative evaluation of performance change for 
this once-per-alrframe-llfe event. The lack of suitable flight dynamic 
load data to conduct other analytical studies of the dynamic load 
dur a % ea H e ^ e1e suggests the need for a flight 

5!? e £ f ! ca 1y designed to obtain propulsion system load 
ill l J?S f L 1 ! ht P^ram suggested could expand the empirical base 
for calibration of the various analytical studies conducted. 



SECTION 2.0 
INTRODUCTION 


The deterioration of a Jet engine's fuel efficiency can be generally 

hJt sSiln' '“.Z" !,0 "^7 >t “ , ' nat1 ? n °( oirfo,,s ’ thermal d 1 start for o? 

J2r 0n P flr ^ 8 '! i flnd ^creases In clearances between rotating blades 
and the corresponding Outer air seals. The h Mi bypass oroDulsfen «!«! 

rf«?Lf 0wo rl n ® toda y' s bedy aircraft fleet' also suffer a shert°term 
R n . P r °cess with losses of about 1 percent In TSFC durlno tho 
VS**'*! ! , f ht s «»t Is undesirable m llaht of current fuel «sts 
high bypass engine deterioration, a study to understand the 
deterioration (Mechanisms has been conducted at Pratt & Whltnev Aircraft 

Ponsorud U p^£* S‘ h °/ s &“** 

r t h , ,." t n l " e consumption of commercial air transports. One'task 
JMn 7 *1^. « s ? en the . utilization of a structural model of the 

loLVon^rfo^nce /et erl o'r a 1 1 oS? POte '’ tia, PffPcts of fl(9ht 1ndu « d 

+lf, V1 i 0U !i ar ) al * tical stud1e s (1) treating flight induced propulsion svs- 

term deter i orat^orf^o'^P a ?hP m f nner su 99 est th ^ primary cause of short- 
J*™ Si ul or i a .!i° n . j be , the increase in operating clearances resulting 

results are n^nilri^^ rot . or/ ^fl interference. The static analysis 
cScIp l In mn jood agreement with engine performance data In the low 
cycle (1 to 300 flight) regime, as shown In Figure 1. 

Although the quasi-steady approach to the flight loads spectrum appears 
to provide an adequate simulation of the clearance Induced detJRSS? 

ture "raising the® D oU 9 imi -? adS F ar * inl r ently time dependent in na- 

possibility of significant dynamic effects in the 

interference process. The transient dynamic Analysis waf 

JTdn 7 / 7 d 7 t fJ nC i re - aSe the , knowledge of the structural behavior of the 
JT9D-7/747 propulsion system and to relate this knowleriop tn coal *»nh 
patterns and performance characteristics exhibited by testing and in 
spection of engines from typical revenue service. 9 

The finite element model of the JT9D-7/747 employed for this ct..Hw wa c 
previously developed by Pratt & Whitney Aircraft fP&WAl and tho Rooi* 

5n°“V AfrP,ane C ° mPany (BCAC) ' Thd dfflmotl th 1s 8 S 

The Intent of the transient dynamic analysis of the JT9D was to examino 

specific load events of a typical flight profile examined in thp onaci 

steady analysis, Figure 3, and determine if the associated dynamic Ifl 

vertirAi !S a l wear . and P erfo ™ance were substantial. Theseevents ar| f a 

t he* Iliac i !t d H USt e . ncodnter and two types of landings. In performing 

P^jectina neak 7 lnad/^th 6y ? mic mp] if Ration factors were used in 
projecting peak loads, therefore, an evaluation of the amplification 


3 


factors was sought. Slnco the dynamic approach requires 
more resources, an evaluation of the relative merits of 
effort * - * oads ffl0( ^ 6 ^ n 9 was a secondary objective of 


considerably 
static versus 
the analysis 


SECTION 3.0 
ANALYTICAL OVERVIEW 


3.1 INTEGRATED PROPULSION SYSTEM MODEL DESCRIPTION 

The NASTRAN finite element model of the OTOD/nocel le structure utilised 
In these studies was Identical to that employed In the steady loads 
analysis. Engine eases are modeled primarily with quadrilateral plate 
elements with enhanced membrane capabilities. Tr Ian- alar plate element*- 
are used In areas of mesh transitioning or Irregular geometry. The 
high-*/ low* pressure rotor spools are beam-lumped mass model duplicates 
of those employed in standard critical speed design analysis of the 

jfgo. 


The model Is substructured by component regions consisting of the 
fan/low-pressure compressor, high-pressure compressor, diffuser/turbine 
cases, the h1gh-/low-pressure rotors, and the inlet cowl, tailcone, and 

wing pylon structures. Since the wing is not modeled, the structure was 

considered symmetric about a vertical plane and one-half of the struc- 

ture is modeled as shown In Figure 4. This half model of the propulsion 
system contains approximately 11,000 degrees of freedom. A summary of 
the model is given in Table I. 

3.2 MODAL ANALYSIS CONSIDERATIONS 

The modal transient approach was selerued as the analysis procedure due 
to its economic superiority Over the r.irect Integration technique for 

this problem. It also provided physical insight into the analysis re- 
garding the frequency spectrum, participation of the various modes in 
the response, and led to the convenient use of modal damping. 

In general the standard NASTRAN transient analysis package was 

utilized. Some modifications were required to the basic NASTRAN com- 
putational procedures, however, due to the following considerations: 

o Desired utilization of structural symmetry 

o Modal damping 

o Time-dependent enforced displacements 

o Nonstandard dynamic reduction procedure 

o Gyroscopic effects 

An overview of the impact of these considerations on the transient 

response formulation follows. 
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• MODAL transient formulation 

15? "Sfj! 4 M anfi l fin i ^ or fu]at1on Is available In many references, (3), 
(4), and wld only be outlined here. The equation to be solved was: ' 

W l»l + CObl * fcl W* M)»b |F(t)| Eqn . i 

where [M] = structure mass matrix 

[C] = structure viscous damping matrix 

[g] = rotor gyroscopic matrix 

DO 5 structure stiffness matrix 

Since the model is substructured the matrices of Equation 1 represent 
enaa^^^ ^ contributions from the various substructures. The zero 
speed, undamped modal matrix of Equation 1 ( 0 N ) was then formed. 


[M] {x} + [K] jx( = {0} 

1*1 = \*o\ e iwt 

+ [k] ) [0 N ] = {of 


Eqn. 2 
Eqn. 3 
Eqn. 4 

I 5 e „'! 0 da j*- matrix ,. *N was m 3Ss normalized and the standard modal 
transformation applied to Equation 1. 


1*1 = [*N]|l| 

t^] T na^N])n(+[» N ] T [c][+ N ] (i| 

+ Ofl] T [G] [♦„] (ill* t* N ]T [K] [* N ] 

* [«N] T )F(t)[ 

['g )■<:+ [«N] T fc][0 N ] |ij(+ [0 N ] t [g] f0 N ]^| 

+ [Vj| >!| • O n ]T |F(t)| 


Eq n . 5 a 


where 


;Kl 
\ 


= generalized mass 
= generalized stiffness 


■M 

f-d 


Eqn. 5b 


Eqn. 6 


Damping in the propulsion system is developed in numerous ways. Bolted 
flanges, case/strut attachments, mount points, design dampers, and In- 
herent structural damping are some examples. In a structure with the 
complexity of the JT9D, an accurate assessment of the discrete dampinq 
matrix [C ] is impractical if not impossible. Consequently, as in many 
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engineering applications, modal damping was employed. The matrix pro- 
duct f 0 w] T [C] fw] thereby Is reduced to diagonal form with 
each entry' proportional to the critical damping factor for the respec- 
tive mode. Modal damping factors may be either obtained from test data 
or assigned values consistent with prior experience. 


3.4 GYROSCOPIC CONSIDERATIONS 


With the Incorporation of modal damping, the equations of motion are 
generally uncoupled, leading to substantial computational savings. The 
gyroscopic behavior of the rotors, characterized by the matrix [GJ , 
produces effects which couple the equations of motion. The physical 
significance of this behavior is that the symmetric and antisymmetric 
components of motion are coupled and the half model used in the quasi- 
steady loads analysis can no longer directly describe the pOwerplant 
structural response. 


A typical stage of the rotor system and the associated analysis coordi- 
nate system is shown in Figure 5. Each stage of the rotor behaves as a 
gyroscope, tending to resist rotation in any direction other than the 
axis of rotation. Mathematically the gyroscopic behavior of each stage 
is defined by Equation 7. 
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Eqn. 7 


where: 


fl = Rotor Rotational Velocity 


To aid in the physical interpretation of the large amounts of displace- 
ment data, it was deemed desirable to maintain the symmetric/antisym- 
metric displacement representation of the structure. This representa- 
tion was preserved by using symmetric and antisymmetric half models of 
the structure and coupling the two representations through the gyro- 
scopic matrix. This technique is applicable only because the coupling 
occurs along the plane of structural symmetry. Using this approach 
Equation 1 may be rewritten as follows: 


'MS 0 ' 


• 1 

*s 


1 

• • 
x A 


r<* 

0 ’ 


*5 

p 

Lo c A „ 


X A 


’K S 

..0 

0 ■ 


*s 

*A 

' ’ 1 


0 


fic« 

0 


i 

*5 


*A 


_eT SA 
f S (t) 

F*(t) 

S and A refer to symmetric and antisvm- 


Eqn. 8 


fE'c dl sp ' oceiient representa t ions V*X SjfTSV' MX 

structural sym- 


sentationf lif? ^ partiti <\ne d into Symmetric and antisymmetric repre- 
sentatTons and the zero speed modal matrix formed, per Equation 4. 
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Eqn. 9 


performed on Tau.TlT « i T a1 ™tn‘ X ’ the modal transformation is 
the JT9D-7 response. ’ ‘ eadin9 to the mathematical representation of 
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Eqn. 10 
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where: 


!!s B symmetric modal participation vector 
^A B antisymmetric modal participation vector 

3.5 DYNAMIC REDUCTION CONSIDERATIONS 


HV1?£ i-TsFz^kWrP'^ 

sta-jssta asr«* s»s£*3s?£ 




necessary to reduce the problem size. The type and level of reduction 
performed can have an adverse effect upon frequency and part 1 cu 1 aly 
modal data. Since the quality of any modal transient analysis Is gov- 
erned by the frequency and modal data used, the dynamic reduction per- 

formed must be done judiciously. 

Several reduction schemes are commonly used, and four types were ex- 
amined for use In this analysis. The types examined were; 

o Mass Lumping followed by Static Condensation 

o Guyan Reduction 

o Distributed Parameter Systems 

o Fixed Boundary Modal Synthesis Techniques 

Both mass lumping and Guyan reduction have been extensively used. Stud- 
ies with these techniques on shell structures have indicated that sub- 
stantial modal errors may be induced even for the low modes of the 

structure when heavy levels of reduction are utilized. Distributed 

parameter systems include the PEER eigenvalue technique (5) and the 
generalized dynamic reduction procedure of MacNeal (7). Both of these 
techniques are relatively new and the latter has shown promising per- 
formance for some engine hardware applications. At the initiation of 
this analysis, however, no evaluation had been made using these tech- 
niques for shell structures. Consequently both techniques were ruled 
out, although subsequent studies have indicated these to be promising 
techniques for future analysis. 


A procedure developed by Craig and Bampton (6) was found to produce 
excellent results for shell structures and to be readily adaptable to 
the NASTRAN program. In this procedure, both physical and modal degrees 
of freedom are used to describe the structure. The modal variables are 
generated by fixing each substructure at its interfaces with other sub- 
structures and extracting the "fixed boundary modes" of the substruc- 
ture. Displacement compatability is then required at the various sub- 
structure boundaries to form the system description, Equation 8. A 
brief theoretical overview of the Craig and Bampton technique follows. 


Each substructure is partitioned into internal and boundary degrees of 
freedom per Equation 11. The boundary set must include all freedoms 
common with other substructures. 
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The fixed boundary frequencies and mode shapes (|<M) of each sub- 
structure are then extracted from equation 12 , 

w? [ M Il] { 0 l} a [kjJ Eqn. 12 

The final displacements of each substructure are then expressed as a 
linear combination of the boundary Induced displacements and the par- 
ticipation of the various fixed boundary modes. 


j X boundary) 

”l 

o" 

Pb 

^ X interna 1 ' 

0 c 

! 0 i 

pi 


where: = substructure boundary displacements 

Pi = participation factors for various substructure fixed 

! ) boundary modes 

0 C >= Displacement of internal degree of freedom due to 
1 boundary displacement, i.e., 

- [Mir 1 • [kib] 

Since each modal degree of freedom ( [Pi] ) represents a physically 
meaningful displacement pattern for the substructure, this technique is 
generally more efficient than Guyan reduction where one freedom typic- 
ally represents a rather arbitrary displacement pattern. Computational- 
ly, the Craig and Bampton technique compares favorably with Guyan re- 
duction. 

In addition to the effectiveness and simplicity of the Craig and Bamp- 
ton technique, the procedure also provides intermediate checks on the 
structural model and analysis. 

o By extracting the fixed boundary modes of the substructures, 

the analyst is given a means of checking the structural model 
early in the assembly phase. 

o The frequency spectrum of the fixed boundary modes provides 

insight into the likely participation of the various modes if 
the harmonic content of the transient loads is known. 

o The fixed boundary mode procedure provides a means for further 
reducing the dynamic problem size. If the fixed boundary fre- 
quencies are considerably higher than the transient forcing 
functions, the fixed boundary modes may simply be restrained 
from participating. 

The need for carefully managing the dynamic reduction procedure is 
illustrated by the simply supported shell of Figure 6 . The results of 
the second bending, first harmonic frequency for this structure are 
shown in Figure 7, comparing the Craig and Bampton procedure with a 
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Guyan reduction scheme similar to the one employed In the original nor- 
mal modes analysis (2). As can he seen, the approach utilized leads to 
a substantial Improvement in the modal data for this case. 

3.6 ENFORCED DYNAMIC DISPLACEMENTS 

In the quasi-steady analysis of the JT9D/747, the wing pylon was at- 
tached to a rigid mount. In dynamic analyses, the powerplant displace- 
ments are due primarily to the elastic response of the aircraft wing. 
Consequently, a procedure must be available whereby time dependent dis- 
placement histories may be enforced on the model.. These displacement 
histories are products of the airframe aeroelastic analysis and are 
defined for all six motions of the wing elastic axis at its Inter- 
section with the pylon centerplane. 

In NASTRAN, enforced dynamic displacements are applied through a "large 
mass" or "acceleration generator". A mass several orders of magnitude 
larger than the structure is assigned to each degree of freedom where 
time dependent displacements are to be applied. The large masses result 
in essentially zer*o frequency, rigid body modes of the structure. These 
rigid body modes may then be excited by forces producing accelerations 
that integrate to the desired displacement histories. 

In the evaluation of this technique prior to this analysis, essentially 
exact results were obtained for simple Spring-mass systems subject to 
the wing aeroelastic displacement histories. However in physical 
modeling of the JT9D/747 and other structures, this method has been 
unable to reproduce the enforced dynamic displacements. The exact error 
involved for the JT9D analysis will be discussed later, however, the 
lack of an accurate technique for enforcing time dependent displace- 
ments lowered the quantitative credibility of the performance 
evaluations obtained for the simulated landing events. 



SECTION 4.0 

DYNAMIC LOADS AND DAMPING DESCRIPTIONS 


The computational costs Involved and limited availability of airframe 
aerolastlc data permitted evaluation of only two flight load conditions 
In the mission profile analyzed In the quasi-steady analysis. The two 
points from the steady loads spectrum which generate the most Interest- 
ing dynamic load conditions are gust encounters and landing touchdowns. 
Both are inherently dynamic In nature and are described In a relatively 
S^rt time period. For these flight conditions, BCAC generated aero- 
e.astlc data for a vertical gust and four different landing sequences. 
The aerolastic data employed reflects the behavior of the wing at the 
outboard engine location. The four landing sequences represent varia- 
tions in aircraft gross weights, fuel loads, and sink rates. With this 
information, three separate dynamic analyses were conducted, 'i he condi- 
tions analyzed corresponded to a once-per-f light vertical gust, a typ- 
ical revenue service landing, and a dnde-per-airplane-life landing. 

4.1 VERTICAL GUST DESCRIPTION 


The vertical gust encounter analyzed was a once-per-f light value of 10 
feet/second at the following flight condition: 


o Altitude 

o Gross Weight 

o Mach Number 

o True Airspeed 


22,100 ft 
730,400 lb 
0.875 

897.8 ft/sec 


Although wing displacement data were available for 1.5 seconds, only 
the first 0.75 second was judged to be valid by the airframe manufac- 
t'^er. For the gust condition, cowl aerodynamic load histories were 
also generated by BCAC and imposed upon the NASTRAN model. 

The coordinate system defining the directions used for wing displace- 
ment histories is shown in Figure 8. The time histories for vertical, 
lateral, and longitudinal displacements along with pitch, yaw, and roll 
values are shown graphically in Figures 9 through 14, respectively. The 
incremental vertical aerodynamic load history created by the chanqinq 
cowl angle of attack is shown in Figure 15. 

4.2 LANDING EVENT DESCRIPTIONS 

Of the four landings described by BCAC, the following two were chosen 
for analysis. 


o Gross Weight 
Fuel Weight 
Sink Rate 


490.000 lb 

43.000 lb 
3 ft/sec 


Gross Weight 
Fuel Weight 
Sink Rate 


449.000 lb 

116.000 lb 
10 ft/sec 


The first ease represents a typical revenue service profile ] art dlng. 
The second a once-per-alrfraitie-llfe event with a hlg.i fuel load and an 
extremely hi gh °s Ink P r ate . This case represents possible the most severe 

dynamic event occurlng during on airframe life. In actual flight tes 

Irtg this condition has been extremely difficult to achieve. 

For the landing events analysed, the Inlet aerodynamic forces _ were 
estimated by 8CAC to be negligible and were neglected. Consequently, 
the dynamic loads for the landing simulations c ft}* J* 
enforced wing displacements. The displacement histories for tne two 
landings are shown in Figures 16 through 27. 

4.3 CALCULATION AND ASSESSMENT OF DAMPING QUANTITIES 

To extract modal damping factors for the propulsion system, accelero- 
meter data from Pratt & Whitney Aircraft's flight testing was ^ VT ^jed. 
Although the flight test aircraft was a B-52, the P^° n ^ engine 
structures for JT9D flight testing are similar to the JT9 J((y A system 
modeled in this analysis. The particular test sequence utilized con- 
Srned verification of the aircraft flight envelope by recording the 
aircraft response to various control surface inputs. Accelerometer data 
from various points on the powerplant were also recorded. 

A hiqh speed transient analysis of the powerplant accelerometer data 
was conducted and the following frequencies identified. 


2.2 

Hz 

25. 

Hz 

5.2 

HZ 

31. 

Hz 

10.2 

HZ 

58. 

HZ 

14.6 

Hz 

60.4 

Hz 


This frequency spectrum is similar to that extracted from the current 
analysis as shown in Table VIII. The trace of the higher modes was of 
insufficient Strength or duration for constructing accurate amplitude 
versus time histories. Consequently a 12 Hz band filter was u ^ z ed 
and an amplitude versus time history was created for the f irst mode of 
the system (2.2 Hz). Using Equation 14 in conjunction with the acceler- 
ometer data, several calculations of damping factors were made and 
averaged. The resulting critical damping factor was 6.06 percent which 
was applied to all modes used in the modal transient analysis. 
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log Aq/Ah 


jl.862 N* + dog A 0 /A n )2|i/Z 


Eqn. 1- 


where: 


6 s critical damping ratio 

A 0 = amplitude of wave at time ti 

A n = amplitude of wave after n th cycle 

N * number of cycles 

Because of the many sources of damping in the propulsion system and 
pylon structures, the value of 6 percent Critical damping was felt to 
be reasonable for the lower modes of the system. In addition, the 
dynamic response Of the propulsion system to the loads being analyzed 
was anticipated to be dominated by the lower modes of the system. Con- 
sequently, an accurate assessment of a typical lower mode damping fac- 
tor was felt to be more important than for higher modes. During the 
coarse of the hard landing analysis, damping studies were performed for 
critical damping factors of 3, 6, and 12 percent with rto substantial 
change in the rub patterns or performance deficit calculated. 


SECTION 6.0 

DYNAMIC RUB STRIP INTERFERENCE MODEL 


8.1 ROTOR/SEAL CLEARANCE CHANGE METHODOLOGY 

HN^rtn r0 thi Sg f by j Wh1 f h ™tor/gag-path seal clearances were calculated 
^!^J?A thi * tr lu ns ana • v "' l§ ,s gh0wr1 schematically In Figure 28. The 
position of the shaft cei r with respect to the case centerline 1$ 
computed for each stage at each time point In the analysis. The change 
In the distance between the deflected and undetected rotor/case loci" 
tlens was then calculated at various points, typically every 12 degrees 

w1trte th L a n a Jn C J r r mf3rence ; Thes t d1stancq changed may ^be compared 
w th the available clearance for each stage at the various cireumferen- 

1 f \ 1 T ! f the clearancGS are opening, closing, or 
' interference (rubbing) has occurred for that time point * The 

available clearance at any time point consists of the P base-line 

* Otv 9 J 6d • 1nt0 t 5 e P. ower P 1ant at a particular performance 
in 1ncu r red P not * t0 the dynamic event as determined 

in tne quasi -steady analysis, and any w^ar incurred as a result of 
previous time points in the dynamic analysis. 

ahhihoU'ii h 16 1n . ter J erence sequence is a complex event, as there is an 

spf 1 U ^nrr.L dyi l° rc& 9ener * ted on the rotor as contact with the 
seal occurs. This force cannot economically be modeled in a linear 

an 2 yS1 !* Consequently, the dynamic interference model is con- 
sidered adequate for small interferences, but not quantitatively 
accurate for large interference values. y 

5.2 TRANSIENT STEP SIZE CONSIDERATIONS 

sldem^oVs" 9 the transient step size there were three governing con- 

o In a modal transient analysis, the step size selected directly 
affects the numerical results through the stability of hiqher 
mode contributions. y 

o The time step must be small enough to track the orbit of each 
rotor in order to provide accurate circumferential clearance/ 
interference information. 

o Economic factors and computing resources. 

f£r „ a ,V Ven t1l " e , 5te J > (At) and the error induced 

l? r (4) Part CiP9t ° 9 m ° de ° f frequenc ^ ( w ) can be expressed by Equation 


Eqrt. 15 


A f 


1 + 


JL 

a 


s * 2 

W At.. 


where: Af * computed amplitude of modal contribution from mode 

with frequency w 

n amplitude of modal contribution from exact solution 
G circular frequency 
a time step size 

It was required that there be less than 6 percent error In a mode with 
a frequency of 100 Hz, and Equation IS was used to yield a time step of 
0.001 second. p 


Af 

A e 

w 

At 


For the dynamic analysis It was anticipated that the propulsion system 
response would be governed by the wing displacements. As can be seen 
from Figures 9 through 14 and 16 through 27, these displacements are 
primarily of relatively low frequency consent when compared to the com- 
puted integrated propulsion system frequency spectrum (Table VIII). 
Consequently, the modes in the 0 to 20 Hz range were expected to have 
* highest participation factors. If, in the modal analysis, modes to 
the 100 Hz level were retained, this would provide a 5 to 1 frequency 
margin over the anticipated primary participation range and would also 
provide adequate coverage for higher modes which might be excited in 
the hard landing, or whose static contributions were required. 


The 0.001 second transient time step provides 8.264 examinations per 
revolution of the high-pressure spool rotor and 17.85 of the low-pres- 
sure spool rotor at cruise conditions. This is considered to be the 
minimum acceptable number of evaluations for the high-pressure spool 
rotor, however, this leads to prohibitly large computer run times. To 
allow this time step to be utilized, velocity and acceleration calcula- 
tions were deleted from the NASTRAN calculations for the various sub- 
structures, saving approximately five hours of computer time per 
dynamic analysis. This procedure limits the amount of data generated 
but does not affect the accuracy of the displacement calculations which 
are used in assessing performance. 
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SECTION 6,0 

JT9P/747 TRANSIENT DYNAMIC RESULTS 


To summarize the Information of the previous sections, the following 
steps were performed in the analysis of the three dynamic flight 
conditions, 3 


o The symmetric and antisymmetric fixed boundary modes of each 
substructure were computed, and the dynamic mass and stiffness 
properties of the various substructures were formed. The rotor 
substructure was reduced using Guyan reduction. 

o Substructure matrlcies were assembled, forr ing Equation 8, and 
the zero speed modes of the Integrated propulsion system were 


0 


0 


0 


Modes of the system up to 100 Hz were retained for a modal 
transient response of the BCAC-supplied dynamic loads with a 
transient step size of 0.001 second. The enforced dynamic dis- 

Efte$l l I l M f1ts °,t t J he wln 9 elastic axis were compared to the 
NASTRAN resultants as a solution check. 


The displacements of the individual substructures were com- 
puted at the various time points and the data used to con- 
struct a performance evaluation via the engine performance 
post-processor computer program. 

An evaluation of the results obtained was conducted, and an 
assessment of the results was made. 


6.1 SUBSTRUCTURE FIXED BOUNDARY MODE SPECTRUMS 

fHJh/hl Xed f b .° 1 ur1c,ar -Y mo . des _ of the cowl, fan, high-pressue compressor, 
turbine, tailcone, and pylon substructures are given in Tables II 
through VII, respectively. As can be seen, the fixed boundary spectrums 
of all substructures exhibit few modes below 100 Hz. 

Frequency and modal convergence studies conducted with the Craiq and 

shown that acceptable propulsion system modal 
convergence is obtained for a fixed boundary frequency cutoff of 
three-to-one over the highest mode of interest. A 100 Hz cutoff value 
for the transient response demands that fixed boundary modes be cal- 
s'/ T a ehi t °fl t \ east , th f . 3 °0 H z level. Examination of Tables II through 
VII shows the lowest fixed boundary cutoff frequency to be 484 Hz for 
the fan symmetric modes. or 


17 


6.? INTEGRATED PROPULSION SYSTEM FREQUENCIES 

oSency of 600 Hz. as were the frequences and mode shapes used In the 
modal transient analysis. 

To investigate the stability of the propulsion svsti cal- 
resoect to various fixed boundary truncation levels, the tr^uencycji 

were* util lied' In P the transient 6 analys™sf°the frequency and *oda, data 
were judged to be numerically acceptable. 

6.3 ENFORCED DYNAMIC DISPLACEMENT BEHAVIOR 

Figures 29 through 34 compare ^ogram^f^^the 

ments representative of this flight event. 

Fioures 35 through 40 compare the input wing displacement histories 

of thedynamlc loading 1s P felt to be accurately analyzed. 

a 4 -v.ai anniuck nf the vertical displacement data for the landing 
simulations Leveed unelplcted high frequency , (above 10 Hz particjpa- 

tion in the aeroelastlc response It .s b iWMl 

suiting 5 from S \he dlfeHn^ .ESS «‘th respect to actual airframe 
response. 

An evaluation of the vertical displacement histories resulting from the 
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loadings on the case and rotors, It Is felt that the conditions 
analyzed are ccr.s c( vdiive with respect to tnose dictated by the aero- 
elastic analysis, 

6.4 VERTICAL GUST RESULTS 


The available clearance for each stage at the time of the gust en- 
counter includes the baseline design clearances at cruise plus any wear 
from Interference which has occurred at previous points In the flight 
prior wear results mainly from aerodynamic loads on the in- 
ana lysTsO) take ’° e ^ climb and are taken from the quasi-steady 


The clrcumferentla 1 rub map resulting from the enee-per-f 1 ight vertical 
gust is shown In Table IX. The resulting TSfC deficits associated with 
this wear are 0.04 percent at cruise and 0.07 percent for the sea level 
take-off point. These 1SFC changes are computed from tables of 
clearance change versus percent change In TSFC (1). 


average 


6.5 REVENUE SERVICE AND HARD LANDING RESULTS 


The circumferential rub patterns resulting from the revenue service and 
hard landing simulations are shown in Tables X and XI, respectively. 
The predicted performance deterioration associated with the revenue 
landing is about 0.2 percent at cruise. The corresponding value for the 
hard landing is about one percent and is believed to be excessive due 
to the lack of a suitable dynamic interference model of the interaction 
between blade tips and outer air seals. 


SECTION 7.0 

ASSESSMENT OF TRANSIENT DYNAMIC ANALYSIS RESULTS 


One problem associated with large scale dynamic analysis Is the genera- 
t i on of a large quantity of data which does not readily yield to a com- 
prehensive evaluation. For this analysis, approximately 99 million 
pieces of data were generated in the process of developing the three 
performance deterioration values for characterizing the gust and 
landing simulations. 

An assessment of the validity of these results can be made through an 
engineering evaluation of various intermediate results of the analysis. 
For the JT9D-7 study, the computed frequency spectrum, performance 
deterioration values and associated seal wear patterns, and the modal 
participation data provides check points for examining the analysis. 

Empirical data for verification of the computed frequency spectrum 
(Table VIII) is somewhat limited. The frequencies obtained, however, 
were in gcOd agreement with the available test data and calculations 
made using the standard Pratt & Whitney Aircraft design analysis tools. 
Table XII "summarizes the available calibration data. 

In addition to these sources of verification, various other computed 
frequencies and mode shapes were found to be in good agreement with 
design studies employing body of revolution, shell of revolution, Or 
beam modeling of miscellaneous engine components. On the basis of the 
available means of correlation, the computed propulsion, system fre- 
quencies were believed accurate, providing confidence in the basic 
elastic and mass modeling of the structures. 

In light of service experience, the once-per-fl ighc vertical gust was 
expected to result in negligible performance loss. The analysis sup- 
ports this experience, predicting a change in TSFC of 0.04 percent from 
the gust simulation. Examination of the predicted seal wear pattern 
from the gust encounter (Table IX) indicates the primary response to.be 
gyroscopic behavior of the fan from the wirg pitching motions. Pitching 
of the low-pressure compressor results in predicted seal wear in the 
last LPC stage. The analysis locations of seal rubs for the fan and 
low-pressure compressor are offset 180 degrees as expected. Inertia 
forces during the gust encounter are not sufficient to induce seal wear 
at the top or bottom of the other air seals, supporting the results of 
the quasi-steady analysis. In summary, the dynamic load induced seal 
wear predicted by the gust simulation is reasonable. 

A performance deterioration of a little over one percent was. calculated 
for the once-per-airplane- 1 ife hard landing event. This flight condi- 
tion is not considered pertinent to the deterioration problem and 
should only be viewed as the upper bound of performance loss which 
could result from extremely rare inertia loadings. The circumferential 
seal wear/rub map for the hard landing (Table XI) indicates rub depths 
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J s b1 9 h a n s . ,D 11s In the fan and 50 mils In the high-pressure 
Hwnam? 6, f ca lu’ 'Interference of this magnitude results In a 

dynamic force on the rotor which Inhibits or may even prevent further 

b1ade . t ,° s@a1 contact. The net result would be considerably less 

seal wear and change In r$FC than predicted by the analysis Althouah 

Important consideration for highly dynamic events such as hard 

an t0 sea ' contTct fore, wa S SU not «'«&? 

auertlv the nr.H(^f»rf a h1J ?. d69 r ee °, f n0,,line8r,t V results. Corise- 
quentiy, the predicted results of seal wear and TSFC change for the 

Suant ItTti Jl 9 rnann^f t10r1 $h ° Uld be v1ewdd oh1y 1n a qual1tat1ve ’ not 

The seal wear/rub patterns resulting from the hard landing again indl- 
cate gyroscopic behavior to be responsible for much of the seal wear. 

the hard . landings are much more complex, however 
as the inertia loads are Sufficient to cause rub at the bottom of the 
case/seals and excite modes of the overhung components such as the 

exDlanatio^ ^ h jJ h "l iressure turbine. A more detailed 

!* P J ? c k f th l se / 1 wear P at tern and rotor behavior for the landina 
events will be made for the revenue landinq. Like the qust encounter 

a slu nd ’ phy a 1° b^seV VS fact ° rS a " d Seal " ear ' a lW to have 


The TSFC debit of about 0.2 percent predicted for the revenue landing 

nation S Sf h lh^ er ct tha t ar1 5 ici P ated » consequently, a more detailed exami- 
" abl $" °L the struct . ural . response to this condition was made. The rub 
map for the revenue landing (Table X) again shows the effects of ovro- 
scopic type (side) motions and rubs. Inertia induced seal wear°is 9 also 
Jnif Ca uL e c d f °I tbe low-pressure compressor and turbine Stages. This re- 
fj 1 * was somewhat surprising as the revenue landing sink rate of 3 

related W seai a wear ht fUel l0ad W ° Uld not be expected t0 produce inertia 


As discussed in Section 6.3, the vertical displacement histories oro- 
Hioh^r NASTRAN enforced dynamic displacement technique result in 

defined bvThe RrAr^I 13 i °f dings J han the dis P 1ad ement histories 
defined by the BCAC aeroelastic analysis. In addition to the in- 

weri S *l<n M r V a *J adin 9J» - pitch rate$ fronrl the aeroelastic analysis 
l° a h 9 ' 6r than anticipated. These two factors increase the gyro- 
scopic behavior of the rotor and provide plausible reasons for the 
revenue service landing seal wear being higher than expected. 

Additional insight into the rotor action and seal wear for the revenue 
landing may be gained by examining the modes of the structure which 
contribute significantly to the structural response! Since the "lode 

16a? e nr(!u! / ripc be a en T? normallzed ' the generalized stiffness (Equation 
tion nrnwiVc atlV t m 1 a i ure °, f the ener 9y of each mode while Equa- 

mode at ^7n,tot“n tta. V * ° f tha a " er9y St0r * d 1n a 9,ven 


[♦] T M [♦]•[«’] Eq '" ,6a 

[ K 

Where Em = medal elastic energy at a given Instant of time. 

The percent of the total elastic energy attributable to each mode of 
the analysis at each time point was tabulated for the revenue landing. 
The modes In which the majority of the system's elastic energy Is con- 
tained are given in Table XIII. The time histories associated with the 
energy involvement of each of these modes are shown in Figure 41 
through 45. 

Although the freqencles of the participating modes are higher than 
anticipated, the sources of excitation for each mode in most cases is 
apparent. The strut lateral bending mode (No. 1) and horizontal fan 
Case rocking mode (No. 16) are typically excited by gyroscopic action 
resulting from wind pitching and the excitation Of modes with fOtdr 
bending. Modes with rotor bending (NoS. 10, 12) are typically excited 
by transient inertia loadings at touchdown. The engine torsion mode 
(No.. 11) is more difficult to rationalize physically but may be excited 
by gyroscopic action in both the vertical and horizontal planes.. Of 
particular interest Is mode 12, which is primarily a low-pressure com- 
pressor case symmetric mode. The excitation of this mode may explain 
the more generalized seal wear occurring in the low-pressure compressor 
for the landing simulation. 


1M.M" 


Eqh. 16b 


[’TO 
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STRUCTURAL MODEL SUMMARY 


Substructure 

Grid Points 

Elements 

Inlet Cowl 

320 

656 

Fan and Low-Pressure Compressor 

359 

604 

High-Pressure Compressor 

464 

516 

Turbine 

587 

826 

Tailcone 

362 

495 

Wing Pylon 

165 

390 

Rotors 

171 

513 

Thrust Frame 

5 

A 



TABLE II 


COWL FIXED BOUNDARY FREQUENCY SPECTRUM 
(Frequencies In Hertz) 


Symmetric 

20.73 

57.60 

114.05 

166.64 

266.14 

304.14 

363.13 

371.13 
456.12 
481.47 

498.00 

507.97 
518.39 
543.23 

560.33 

571.33 
574.55 

597.97 
613.85 

634.14 

638.15 
639.11 
656.27 
696.37 
728.66 
764.69 
864.42 
871.61 

911.26 
935.22 
970.45 
987.18 

1094.00 

1164.26 

1206.26 
1238.74 
1302.46 
1330.65 
1360.57 


Antisymmetric 

14.75 

369.51 
275.84 
288.05 
372.62 
423.27 
450.76 
468.94 
496.96 
511.44 

527.57 

543.25 

555.58 

569.16 

622.51 

626.16 
670.36 
691.30 
783.50 
786.80 

828.65 
873.47 
903.05 

1020.83 

1037.65 
1157.57 
1175.59 

1195.84 
1285.32 
1295.64 
1301.69 
1336.78 
1343.82 

1366.25 
1545.02 
1565.99 
1597.67 
1640.61 
1727.54 


1436.84 

1482.10 

1494.50 

1550.27 

1560.19 

1635.95 

1649.09 

1758.25 
1910.88 
2018.04 
2040.94 
2128.87 
2148.81 

2184.35 

2270.35 
2429.67 

2632.26 


1741.61 

1776.02 

1915.32 

2021.95 

2129.97 

2143 . 02 . 

2185.00 
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TABLE IU 


PAN FIXED BOUNDARY FREQUENCY SPECTRUM 
(Frequencies In Hertz) 


Symmetric 


Antisymmetric 


95.05 

330.68 

93.88 

352.36 

103.15 

332.89 

100.49 

355.48 

126.73 

338.26 

126.77 

357.05 

134.17 

353.39 

137.61 

358.58 

161.13 

353.99 

159.75 

369.99 

175.98 

358.98 

176.12 

372.08 

176.44 

362.18 

177.50 

386.34 

178.48 

369.97 

178.33 

402.78 

183.08 

375.59 

181.68 

404.00 

185.77 

393.73 

184.12 

407.22 

186.40 

402.79 

ia&.b2 

411.34 

187.77 

404.02 

187.97 

423.87 

189.52 

407.55 

189.69 

430.99 

189.91 

407.74 

190.83 

436.36 

191.46 

409.74 

204.03 

438.47 

204.60 

414.32 

204.94 

439.92 

204.92 

420.05 

205.13 

447.96 

205.92 

429.89 

208.17 

453.66 

208.80 

430.32 

233.46 

462.57 

220.78 

434.85 

241.68 

462.95 

233.47 

436.38 

241.95 

468.24 

241.83 

439.21 

244.01 

484.34 

242.61 

439.93 

245.59 

485.79 

244.96 

445.65 

246.71 

494.66 

247.22 

448.20 

263.19 

498.00 

249.14 

456.45 

267.78 

502.42 

265.53 

461.81 

271.40 

506.83 

268.36 

462.81 

273.62 

523.56 

272.88 

463.19 

275.92 

530.78 

273.25 

464.62 

276.71 

532.57 

276.53 

471.38 

278.97 

537.79 

277.56 

480.59 

280.39 

539.86 

280.88 

484.53 

289.19 

541.13 

281.34 


289.29 


281.78 


294.19 


289.14 


301.37 


290.05 


311.63 


294.17 


316.82 


309.18 


330.70 


313.94 


332.15 


316.03 


337.34 


320.31 


338.69 
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TABLE IV 


HPC FIXED BOUNDARY FREQUENCY SPECTRUM 
(Frequencies In Hertz) 


Antisymmetric - 


126.70 

669.95 

157.09 

712.67 

163.60 

671.59 

174.38 

730.60 

194.04 

682.51 

216.02 

737.95 

214.47 

699.42 

234.63 

745.09 

250.76 

703.77 

253.45 

750.58 

282.73 

720.45 

307.98 

752.61 

300.42 


314.16 


306.59 


335.83 


331.18 


347.58 


336.41 


384.96 


360.08 


388.35 


373.57 


411.36 
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TABLE V 


TURBINE FIXED BOUNDARY FREQUENCY SPECTRUM 
(Frequencies In hertz) 

Symmetric Antisymmetric 


as . a ? 

457.10 

616.96 

43.69 

464.57 

677.55 

47.96 

463.79 

621.02 

77.11 

473 . 9 ? 

677.99 

76.12 

465.34 

640.82 

88.13 

477.92 

679.89 

86*64 

477.02 

646.20 

98.10 

485.10 

684.26 

107.66 

480.14 

648.08 

108.58 

487.20 

694.46 

119.29 

483.89 

653.23 

120.25 

490.60 

703.39 

126.68 

484.36 

669.54 

128.80 

491.88 

709.40 

128.99 

488.83 

672.32 

128.99 

492.90 

713.51 

129.38 

490.73 

676.79 

129.05 

499.18 

739.60 

133.55 

491.69 

677.92 

134.17 

499.81 

743.45 

134.51 

493.72 

679.69 

134.93 

505.75 

744.17 

140.22 

494.37 

680.42 

141.06 

508.27 

751.93 

142.91 

499.27 

697.16 

143.33 

512.84 

762.78 

143.34 

501.11 

709.71 

150.28 

516.17 

775.90 

153.71 

508.22 


153.90 

521.60 


162.36 

510.68 


162.55 

523.27 


165.85 

512.02 


195.63 

528.24 


192.31 

516.05 


216.21 

529.24 


207.90 

520.32 


223.57 

535.85 


215.75 

523.60 


231.99 

545.42 


223.67 

528.60 


243.06 

545.98 


242.76 

535.91 


250.24 

548.33 


249.96 

540.66 


256.95 

550.62 


257.07 

543.37 


267.72 

560.02 


264.72 

545.63 


278.22 

568.85 


267.88 

546.07 


282.62 

573.25 


277.17 

548.99 


293.41 

577.88 


284.01 

550.95 


305.52 

581.08 


301.61 

552.28 


309.12 

586.21 


306.29 

564.53 


316.02 

593.25 


309.68 

569.34 


350.19 

593.47 


349.90 

573.33 


350.36 

593.73 


361.15 

578.06 


362.60 

595.95 


363.55 

592.01 


363.80 

598.92 


364.70 

593.53 


365.67 

600.38 


367.85 

594.54 


366.43 

606.05 


386.56 

598.74 


389.84 

606.46 


388.47 

599.31 


390.81 

610.26 


391.75 

600.99 


418.62 

642.34 


401.74 

601.29 


426.50 

644.09 


417.93 

605.14 


428.87 

659.91 


428.77 

605.67 


454.85 

669.31 


432.25 

607.34 


462.70 

674.67 
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TABLE VI 


TAILCONE -FIXED BOUNDARY FREQUENCY SPECTRUM 
(Frequencies In Hertz) 


Symmetric 

85.12 
222. 55 
329.67 
367.02 
373.20 

392.37 

569.44 

668.15 
721.65 

882.33 
907.50 

951.22 

990.76 
1099.93 
1152.42 

1267.76 

1279.99 

1317.25 
1337.23 
1355.81 

1377.19 
1404.96 

1436.15 
1451.07 
1471.91 

1488.16 

1496.44 
1526.41 

1539.25 

1625.22 
1783.56 
1921.58 

2023.64 

2136.77 
2202.30 

2236.37 

2402.99 

2466.19 
2531.69 

2786.60 
2828.89 

2873.65 


Antisymmetric 

85.07 
226.68 

328.61 

336.62 
369.78 

376.43 
575.93 

687.39 

721.43 
917.70 
985.36 

100.33 

1100.39 

1158.61 

1192.38 

1277.04 
1285.96 
1324.47 

1346.05 
1362.88 
1379.53 

1418.33 

1461.22 
1467.13 

1477.60 
1502.00 
1513.83 
1529.02 

1920.27 

1930.52 

1975.27 
2049.29 
2196.42 
2248.12 

2272.99 

2403.52 
2878.56 

2942.61 

2947.40 

3126.25 

3304.45 

3433.23 


2886.15 

3186.83 
3275.93 
3436.48 

3492.84 
3933.30 
4259.05 
5322.21 


3612.79 

3932.11 

4316.98 

4413.85 

4707.55 

4943.55 
5424.05 
6576.91 





Mode 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


table VIII 

integrated propulsion system frequencies 

* Frequency Test Results 
S or A (Hz) (Hz) _ Description 


A 

2.5190 

2.4 

S 

4.4986 

4.6 

A 

6.8272 

6.5 

S 

11.115 


A 

12.836 


A 

14.497 


S 

14.816 


A 

19.123 


S 

20.089 


S 

24.510 


A 

27.489 


S 

29.024 


s 

31.809 


A 

32.499 


S 

35.126 


A 

37.833 


S 

45.539 


A 

45.556 


S 

46.648 


A 

51.473 


S 

53.351 


A 

54.059 


S 

56.636 

56. 

S 

62.207 


A 

63.216 


S 

70.002 


S 

70.896 


A 

72.600 


S 

74.723 


A 

75.101 


A 

75.865 


A 

76.513 


S 

77.930 


A 

83.489 


s 

86.212 


S 

87.184 


A 

87.775 


A 

88.225 


S 

89.475 


A 

90.955 



First Pylon Lateral 
First Pylon Vertical 
Strut Torsion 
Thrust (Y) Frame, Engine 
Axial 

Engine Roll 

Fan Case Rocking - Side 
Fan Case Rocking - Vertical 

Fan Thrust Reversers 
Higher Pylon Bending With 
Low Rotor Motion 
Case Torsion, Low Rotor 
Bending 


Core Engine Bending 


Tailcone Wagging, Low Rotor 
Bending 


★ 
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S = symmetric mode 
A = antisymmetric mode 


TABLE VIII (Cont'd) 

INTEGRATED PROPULSION SYSTEM FREQUENCIES 


Mode 

* 

S or A 

Frequency 

(Ha) 

41 

S 

93.189 

42 

A 

93.732 

43 

S 

94.828 

44 

A 

95.689 

45 

A 

97.376 

46 

A 

98.838 

47 

S 

99.287 

48 

A 

101.64 

48 

S 

102.75 

SO 

A 

103.75 

51 

S 

104.26 

52 

A 

106.65 

53 

5 

107.65 

54 

S 

107.79 

55 

A 

108.57 

56 

$ 

114.22 

57 

A 

115.61 

58 

S 

116.06 

59 

S 

116.03 

60 

A 

118.14 

61 

S 

119.31 

62 

A 

120.27 

63 

A 

120.35 

64 

A 

122.94 

65 

S 

125.73 

66 

S 

126.28 

67 

A 

126.77 

68 

S 

126.80 

69 

s 

126.90 

70 

s 

127.67 

71 

A 

128.80 

72 

A 

128.98 

73 

A 

128.99 

74 

S 

129.02 

75 

A 

129.04 

76 

S 

129.41 

77 

A 

130.71 

78 

A 

132.54 

79 

S 

133.57 

80 

A 

134.17 


Test Results 



Description 


S - symmetric mode 
A = antisymmetric mode 
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TABLE IX 



i 


V 


CIRCUMFERENTIAL RUB PATTERN FOR 
QNCE-PER-FLIGHT VERTICAL GUST 
(Rub Depths In mils) 


Engine Stages 


0 P 8 

6 A 

9 N L 

» P 

6 C 

6 

8 * 

I I I-— 

0. 0. o. 

U. 0. o. 

24. 0. 6, 

36. 3. 0. 

46. 9, o. 

66. 15. 6, 

72. ZB. o. 

86.. sa, o. 

66. 65. 6. 

106. 3S . 0. 

160. 19. 0. 

1 32. 6. 6. 

146. 0. 0. 

1S6. 0. 0. 

168. 0. 0. 

160 . 0 . 0 . 

192. 0. 0. 

206. 0. 0. 

216. 0. 0. 

226 . 0 . 0 . 

240. 0. 0. 

252. 0. 0. 

266. 7. o. 

276. 16. 0. 

268. 6. 0. 

300. 3. 0. 

312. 1. 0. 

326. 0. 0. 

336. 0. 0. 

346. 0. 0. 



♦Circumferential location In degrees, viewed (clockwise) from front of 
the engine; 0 degrees Is at top of the engine. 




TABLE XI 




CIRCUMFERENTIAL RUB PATTERN FOR 
HARD LANDING SIMULATION 
(Rub Depths In mils) 


Engine Stages 



OH 

64 

: 


I t - 

1- 

— I- 

— I- 

— 

— 1- 

X- 

— I- 

X- 

. 1 

i 

.... 

-I — 

-I — 

-X I- 

X- 

»«■«• 

6. 

0. 

0. 

0. 

0. 

0. 

6. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

. o. 

0. 

0. 

0, 

0. 

0. 

0. 

12. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

24. 

0. 

0. 

0. 

0, 

0. 

0. 

0. 

0. 

0. 

b. 

0. 

0. 

0. 

0. 

0. 

1. 

0. 

0. 

0. 

0. 

36. 

0. 

0. 

0. 

12. 

0. 

0. 

0. 

0. 

6. 

0. 

0. 

0. 

6. 

0. 

0. 

15. 

1. 

0. 

0. 

0. 

40. 

0. 

0. 

0. 

36, 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

5. 

27. 

14. 

0. 

0. 

0. 

60. 

0. 

0. 

1. 

69. 

0. 

0. 

0. 

0. 

«. 

0. 

0. 

0. 

0. 

o. 

15. 

39. 

31. 

7. 

0. 

0. 

72. 

16. 

0. 

10. 

66. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0, 

22. 

46. 

53. 

21. 

0. 

0. 

64. 

72. 

0. 

19. 

68. 

0. 

0. 

0. 

0. 

h 

0. 

b. 

0. 

0. 

0. 

21 . 

48. 

67. 

28. 

b. 

0. 

96. 

136. 

4. 

21. 

78. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

b. 


0 

15. 

42. 

67. 

29. 

0. 

0. 

106. 

165. 

14. 

17. 

61 . 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

b. 

6. 

31. 

55. 

22. 

0. 

b. 

120. 

161. 

30. 

IS. 

45. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

15. 

37. 

6. 

0. 

0. 

132. 

136. 

41. 

27. 

39. 

0. 

0. 

0. 

0. 

0. 

0. 

b. 

0. 

b. 

o. 

0. 

0. 

19. 

0. 

0. 

0. 

144. 

77. 

46. 

45. 

42. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

6. 

0. 

0. 

5. 

0. 

0. 

0. 

156. 

27. 

43. 

58. 

50. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

. o. 

0. 

0. 

0. 

0. 

0. 

0. 

160. 

0. 

35. 

64. 

56. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

6. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

166. 

6. 

20. 

63. 

60. 

0. 

6. 

0. 

0. 

0. 

6. 

0. 

0. 

0. 

b. 

b. 

0. 

0. 

0. 

0. 

0. 

192. 

0. 

37. 

62. 

57. 

0. 

0 

0. 

0. 

0. 

0. 

0. 

0. 

b. 

0. 

6. 

0. 

0. 

0. 

0. 

0. 

204. 

19. 

47. 

53. 

47. 

0. 

0. 

0. 

0. 

0. 

b. 

0. 

0. 

0. 

. o. 

0. 

0. 

0. 

0. 

0. 

0. 

216. 

49. 

50. 

35. 

36. 

0. 

6. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

. b. 

0. 

0. 

12. 

0. 

b. 

0. 

226. 

77. 

47. 

14. 

33. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

6. 

0. 

11. 

31. 

5. 

0. 

0. 

240. 

92. 

37. 

0. 

32. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

14. 

26. 

51. 

22. 

b. 

o. 

252. 

60. 

16. 

5. 

45. 

6. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

26. 

39. 

70. 

36. 

0. 

0. 

mi . 

43. 

4. 

11. 

70. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

o. 

31. 

47. 

80. 

44. 

6. 

0. 

276. 

13. 

0. 

17. 

67. 

0, 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

. o. 

29. 

5b. 

80. 

44. 

0. 

0. 

266. 

0. 

0. 

19. 

66. 

0. 

0. 

0. 

0. 

o. 

0. 

0. 

0. 

0. 

o. 

22. 

49. 

66. 

36. 

0. 

0. 

300. 

0. 

0. 

16. 

74. 

0. 

0. 

0. 

0. 

0. 

b. 

b. 

o. 

b. 

0. 

11. 

35. 

49. 

20. 

0. 

0. 

312. 

0. 

0. 

7. 

53. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

19. 

30. 

3. 

0. 

0. 

324. 

0. 

0. 

2. 

19. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

o. 

0. 

0. 

3. 

15. 

0. 

0. 

0. 

336. 

6. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

346. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

b. 

b. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 


*Cireumferential location In degrees, viewed (clockwise) from front of 
the engine; 0 degrees is at top of the engine. 
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TABLE. XII 


CALIBRATION OF COMPUTED FREQUENCY SPECTRUM 



NASTRAN Analysis 
(Hertz) 

Test(T) or Analysis(A) 
(Hertz) 

Strut Lateral Bending 

2.52 

2.40(T) 

Strut Vertical Bending 

4.50 

4.60(T) 

Strut Torsion 

6.50 

6.80(T) 

Tailcone Vertical Wagging 

56.60 

56.00(T) 

Fan Case Rocking/Low Rotor 

24.51 

22.83(A) 

LPC Case Bending/Low Rotor 

2.9.02 

30.17(A) 
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TABLE XIII 

HIGH ENERGY MODES FOR TYPICAL REVENUE SERVICE LANDING 


Mode No. 

Frequency 

(Hertz) 

General Description 

1 

2.51 

Strut Lateral Bending 

10 

24.51 

Vertical Fan Rocking, Low-Pressure Spool 
Rotor Bending 

11 

27.49 

Engine Torsion, Low-Pressure Spool Rotor 
Bending 

12 

29.02 

Low-Pressure Compressor Case Vertical 
Bending, Low-Pressure Spool Rotor Bending 

16 

37.833 

Core Engine Bending, Lateral Fan Case 
Rocking 
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TT 


■5W 


T 






° JTdD-7A 

♦5.0 - n JT9D-7A ISP) 
a JT9D-7F 
OJT9D-20 

% ATSFC ,4 '° 
relative ♦3.0- a nastran 
to new o 


O “ 

♦1.0- 

<S O 


50 100 150 200 

Flight cycles 


250 300 


Fl9Ure 1 A°J? JsifSf sw!* < m' -' l PecformBnee Deterioration Based on 

Analysis of Steady Flight Loads with Fleet Performance Data 



Figure 2 JT9D-7/747 Integrated NASTRAN Finite Element Structural Adel 
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Ground runs T.O. Cruise 


Figure 3 


Speed 

checks 



Stall 
checks. 

Descent-' Approach 


\ 


Landing 


Typical Flight Profile Analyzed in the Steady Loads Analysis 



re 4 JT9D-7 Finite Element Model Displaying Use of Structural 
Symmetry 







RADIAL DEFLECTION (RELATIVE* 



figure 7 


M6de Shape Comparisons for Simply Supported Cylindrical 
Shell of Figure 6 



Figure 8 Wing Displacement Coordinate System 
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i iriE ^ SEC ) 

Figure 9 Wing Vertical Displacement History for Once-Per-Flight 
Vertical Gust 



tike ^sec) 


Fiqure 10 Wing Lateral Displacement History for Once-Per-Flight 
Vertical Gust 













Figure 16 Wi rig Vertical Displacement History for Typical Revenue 
Service Landing 
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Figure 17 Wing Lateral Displacement History for Typical Revenue 
Service Landing 



Figufe 18 



Wirtg Axial Displacement History for Typical Revenue 
Service Landing 
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Figure 19 Wing Pitch History for Typical Revenue Service Landina 
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Figure 20 Wing Yaw History for Typical Revenue Service Landing 





figure 22 


Figure 23 



Wing Vertical Displacement History for Hard Landing 



Wing Lateral Displacement History for Hard Landing 
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Figure 24 Wing Axial Displacement History for Hard Landing 
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Figure 25 Wing Pitch History for Hard Landing 
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Figure 26 Wing Yaw History for Hard Landing 
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A. ROTOR/CASE UNDEFLECTED POSITION AT TIME = t Q 


CLEARANCE 



TYPICAL CLEARANCE CHANGE 

CALCULATION AT EACH RADIAL 

LOCATION (r • r = 5 RUB) 
c r 


DIRECTION OF ROTOR AND 
CASE DEFLECTION 


6. ROTOR/CASE DEFLECTED POSITION AT TIME = t 0 +At 



ROTOR 


AVAILABLE CLEARANCE INCREASED 
BY PREVIOUS RUB (NON UNIFORM 
CIRCUMFERENTIALLY) 


C. ROTOR/CASE DEFLECTED POSITION AT SUBSEQUENT TIME POINT 


Figure 28 


Dynamic Clearance Change Calculation 


Procedure 
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V. 




wlilfl meplflcemertt Input 
«■. «— NASTftAN progrftrri Output 



Figure 29 Wing Vertical Displacement Evaluation for Vertical Gust 


SCAC wing displacement Input 

NASTfiAN program Output 



RFPRODUCIBILITY OF Till 

SSm. PM5B IS ««* 


Figure 30 Wing Lateral Displacement Evaluation for Vertical Gust 
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■g r-wrsnr 


Figure 31 


— ’ — BCAC wing dlBpIflcemqnt Mipwt 
— — MA8TRAN prn^nm output 



o*.go o!ao 
T I ht‘( SEC ) 


Wing Axial Displacement Evaluation for Vertical Gust 


QCAC wing displacement Input 

NASTRAN program output 



Figure 32 Wing Pitch Displacement Evaluation for Vertical Gust 



— ^—BCAC wlrifl dlafllAGflmortt Ifipwt 
— ^MASTftAN program output 
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figure 33 Wing Yaw Displacement Evaluation for Vertical Gust 


&CAC wing displacement Input 

— —NASTRAN program output 



Figure 34 Wing Roll Displacement Evaluation for Vertical Gust 





— fiCAC wing displauamuni input 
MAST RAM program output 


Revenue Service Uncling .. Hard Landing 



Figure 35 Wing Vertical Displacement Evaluation for Landing Events 


— — BCAC wing displacement Input 
■“ — NASTRAN program output 

Revenue Service Landing Hard Landing 



Figure 36 Wing Lateral Displacement Evaluation for Landing Events 



— -,r-6CAC wing dlaplacemailt input 

Revenue Service Landing nabthan prmwm pm*™ Hard Land i n g 



Figure 37 Wing Axial Displacement Evaluation for Landing Events 


— HCAC wing displacement Input 

Revenue Service Landing --NASTrtANpr O0 , flm#o , PU , Hard Land1r1g 



Figure 38 Wing Pitch Displacement Evaluation for Landing Events 






BCAC wing displacement input 

Revenue Service Landing nabtran pmarnm output Hard Landing 



Figure 39 Wing Yaw Displacement Evaluation for Landing Events 


— ^ ~ BC AC wing dlsplacumum Input 

Revenue Service Landing * —-nastran ptogmm output Hard Landinn 



Figure 40 Wing Roll Displacement Evaluation for Landing Events 
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il 



Figure 41 Energy Level Associated with Mode 1 for the Revenue 
Service Landing 


MOOt r 



'V OK j'fi r 


Figure 42 Energy Level Associated with Mode 10 for the Revenue 
Service Landing 
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Figure 43 Energy Level Associated with Mode 11 for the 
Service Landing 
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Figure 44 Enr.-gy Level Associated with Mode 12 for the 
Service Landing 
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Figure 45 
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Energy Level Associated with Mode 16 for the Revenue 
Service Landing 
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APPENDIX 
LIST OF SYMBOLS 


A 

m 

Amplitude of mode 

E 

R 

Elastic energy of system 

g 

s 

Gravity 

i 

s 

Mass moment of inertia of stage 

M 

e 

Moment 

N 

= 

Number of Cycles 

h 

s 

Substructure boundary displacements 

pi 

s 

Participation factors for various 
boundary modes 

TSFC 

r 

Thrust specific fuel consumption 

At 

= 

Incremental time step 

0 

r 

Rotation 

n 

= 

Rotor rotational velocity 


:r 

Critical damping ratio 

n 

f 

Displacement of internal degree 
boundary displacement, i.e., |0 C ; 

O) 

= 

Circular frequency 

[c] 

= 

Structure viscous damping matrix 

[ G ] 

s 

Rotor gyroscopic matrix 

[«] 


Structure stiffness matrix 

!> 

"i 

.] . 

Generalized stiffness = [^2 N J 

[ M ] 


Structure mass matrix 


] ■ 

Generalized mass = 

[♦] 

= 

Zero speed modal matrix 

H 

s 

Displacement vector (in) 

{F(» 

>} ■ 

Forcing function (lbs) 

W 

a 

Modal participation vector 
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Subscripts 


A 

B 

Anti "Symmetric 

b 

P3 

Boundary 

e 

* 

Exact 

f 

e 

Frequency 

I 

s 

Internal 

m 

ft 

Modal 

N 

= 

Normalized 

rt 

= 

Ntb cycle 

0 

s 

Initial 

P 

= 

Polar 

S 

s 

Symmetric 

t 

= 

Time step size 

x,y,z 

= 

Coordinate direction 

Matrix 

Symbols 


[] T 

s 

Matrix transpose 

[ ] _1 

s 

Matrix inverse 

• • 

[ ] 

= 

Acceleration 

[*] 

= 

Velocity 

C^] 

= 

Diagonal matrix 


tt'IPRODUClBlUTY OF '111b 
ORIGINAL PACK IB P0 ( Hl 


65 






W- 


W: 








"■ «>";«? ‘“' 


rSF7 


JV 


